Yang— Baxterizat ions, Universal Quantum 
Gates and Hamiltonians 



Yong Zhang 1 

Institute of Theoretical Physics, Chinese Academy of Sciences 
P. 0. Box 2735, Beijing 100080, P. R. China 

Louis H. Kauffman 2 
Department of Mathematics, Statistics and Computer Science 
University of Illinois at Chicago 
851 South Morgan Street, Chicago, IL, 60607-7045 

Mo-Lin Ge 3 

Nankai Institute of Mathematics, Nankai University 
Tianjin 300071, P. R. China 



Abstract 

The unitary braiding operators describing topological entanglements can be 
viewed as universal quantum gates for quantum computation. With the help 
of the Brylinski's theorem, the unitary solutions of the quantum Yang-Baxter 
equation can be also related to universal quantum gates. This paper derives uni- 
tary solutions of the quantum Yang-Baxter equation via Yang-Baxterization 
from the solutions of the braid relation. We study Yang-Baxterizations of 
the non-standard and standard representations of the six-vertex model and 
the complete solutions of the non-vanishing eight-vertex model. We construct 
Hamiltonians responsible for the time-evolution of the unitary braiding opera- 
tors which lead to the Schrodinger equations. 
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1 Introduction 



There are natural relationships between quantum entanglement and topolog- 
ical entanglement :2.. Topology studies global relationships in spaces, and how 
one space can be placed within another, such as knotting and linking of curves in 
three-dimensional space. One way to study topological entanglement and quan- 
tum entanglement is to try making direct correspondences between patterns of 
topological linking and entangled quantum states. One approach of this kind 
was initiated by Aravind [Sj, suggesting that observation of a link would be mod- 
elled by deleting one component of the link. But this correspondence property of 
quantum states and topological links is not basis independent 

A deeper method (we believe) is to consider unitary gates R that are both 
universal for quantum computation and are also solutions to the condition for 
topological braiding. Such ^-matrices are unitary solutions to the Yang-Baxter 
equation (the braid relation). We are then in a position to compare the topologi- 
cal and quantum properties of these transformations. In this way, we can explore 
the apparently complex relationship among topological entanglement, quantum 
entanglement, and quantum computational universality. This way has been ex- 
plored in a series of papers @1 |Hl El O IH1 El [101 [TT1 H^J . 

The present paper derives unitary solutions of the Quantum Yang-Baxter 
Equation (QYBE) via Yang-Baxterization and explores the corresponding dy- 
namical evolution of quantum entanglement states. The solutions to the QYBE 
that we derive by Yang-Baxterization contain a spectral parameter x, and hence 
do not, except in special cases, give representations of the Artin braid group. 
These new solutions are unitary, and they do give useful quantum gates. Thus 
we show in this paper that the full physical subject of solutions to the quan- 
tum Yang-Baxter equation (including the spectral parameter x) is of interest for 
quantum computing and quantum information theory. 

The plan of the paper is organized as follows. In the second section, a unitary 
braiding operator is regarded as a quantum entanglement operator and further 
a universal quantum gate. To describe the dynamical evolution of the unitary 
braiding operator R, Yang-Baxterization is used to solve the QYBE. In the third 
section, the unitary solutions of the QYBE for the non-standard and standard 
representations of the six-vertex model are obtained via Yang-Baxterization. In 
the fourth section, the complete unitary solutions of the QYBE for the non- 
vanishing eight-vertex model are obtained via Yang-Baxterization. In the fifth 
section, with the Brylinksi's theorem JH], all unitary i?(x)-matrices presented in 
this paper are recognized to be universal quantum gates for most x-values. In the 
sixth section, the Hamiltonian determining the time evolution of quantum state is 
constructed with the unitary ^(x)-matrix. In the seventh section, as an example, 
the CNOT gate is constructed in terms of the unitary ^-matrix (_R(a;)-matrix) 
and local unitary transformations. To conclude, remarks on our work are made. 
In Appendix A, a pragmatic introduction to Yang-Baxterization is presented. 
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2 The QYBE in quantum entanglements 



This section presents basic elements underlying our work. The Braid Group Rep- 
resentation (BGR) and the QYBE are described in the sense of studying quantum 
entanglements. The Brylinski's theorem ^3] plays the key role in relating a uni- 
tary i?-matrix (i?(x)-matrix) to a universal quantum gate. Yang-Baxterization 
derives the corresponding .R(x)-matrix from the ^-matrix and makes it possible 
constructing the Hamiltonian determining the evolution of the unitary braiding 
operator (i?-matrix). 

2.1 The BGR in quantum entanglements 

Braids are patterns of entangled strings. A braid has the form of a collection of 
strings extending from one set of points to another, with a constant number of 
points in each cross section. Braids start in one row of points and end in another. 
As a result, one can multiply two braids to form a third braid by attaching the 
end points of the first braid to the initial points of the second braid. Up to 
topological equivalence, this multiplication gives rise to a group, the Artin braid 
group B n on n strands, which is generated by < i < n — 1}. 

The group B n consists of all words of the form b^b^.-.b^ modulo the rela- 
tions: 

bib i+ ibi = bi+ibibi+i, 1 < i < n - 1, 

bibj = bjbi, \i - j\ > 1. (1) 

Each braid is, in itself, a pattern of entanglement. Each braid is also an oper- 
ator that operates on other patterns of entanglement (braids) to produce new 
entanglements (braids again). 

The analogy between topological entanglement and quantum entanglement 
(from the point of view of braids) means the association of a unitary operator with 
a braid that respects the topological structure of the braid and allows exploration 
of the entanglement properties of the operator. In other words, we propose to 
study the analogy between topological entanglement and quantum entanglement 
by looking at unitary representations of the Artin braid group. The main point 
for the exploration of the analogy is that (from the point of view of a BGR) each 
braid is seen as an operator rather than a state. See Fig. 1. 

Consider representations of the Artin braid group constructed in the following 
manner. To an elementary two strand braid there is associated an operator 

R:V®V — ► V ® V. (2) 

Here V is a complex vector space, and for our purposes, V will be two dimensional 
so that V can hold a single qubit of information. One should think of the two 
input and two output lines from the braid as representing this map of tensor 
products. Thus the left endpoints of R as shown in Fig. 1 and Fig. 2 represent 
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the tensor product V ® V that forms the domain of R and the right endpoints of 
the diagram for R represent V ® V as the range of the mapping. In the diagrams 
with three lines shown in Fig. 2, we have mappings from V ® V ® V to itself. 

The identity shown in Fig. 2 is called the Yang-Baxter Equation (the braid 
relation), and it reads algebraically as follows, where Id denotes the identity 
transformation on V. 

(R®Id)(Id®R)(R® Id) = (Id®R){R® Id) (Id ® R) . (3) 

This equation expresses the fundamental topological relation in the Artin braid 
group, and is the main requirement for producing a representation of the braid 
group by this method. 

2.2 The Brylinski's theorem and ^-matrix 

A two-qubit gate G is a unitary linear mapping G : V ® V — ► V (8> V where V 
is a two complex dimensional vector space. A gate G is said to be entangling if 
there is a vector 

\a(3) = \a) ® \/3) € V (8) V 

such that G\a@) is not decomposable as a tensor product of two qubits. Under 
these circumstances, one says that G\a0) is entangled. 

In ! the Brylinskis give a general criterion of G to be universal (in the 
presence of local unitary transformations). They prove that a two-qubit gate G 
is universal if and only if it is entangling. 

Here is the specific i?-matrix that we shall examine 



R 



(a 
















d 








c 








V o 








b J 



(4) 



Id® Id® Id 





R® Id 



Figure 1: A braiding operator .R<8> Id. 
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{R® Id)(Id® R)(R® Id) = (Id®R)(R® Id) (Id <g> R) 
Figure 2: The Yang-Baxter equation (the braid relation). 



where a, b, c, d can be any scalars on the unit circle in the complex plane. Then 
R is a unitary matrix and it is a solution to the Yang-Baxter equation 

The point of this case study is that R, being unitary, can be considered as 
a universal quantum gate and since R is the key ingredient in a unitary repre- 
sentation of the braid group, it can be considered as an operator that performs 
topological entanglement. We shall see that it can also perform quantum entan- 
glement in its action on quantum states. The ^-matrix can also be used to make 
an invariant of knots and links that is sensitive to linking numbers [2]. 



\<t>) 



£(M®M) = |0) 



Figure 3: The braiding operator R as an entangling operator. 



Consider the action of the unitary transformation R on quantum states. We 
have 

£100} = o|00), £|01) = c|10), 

R\W) = d\01), R\ll) = b\U). (5) 

Here is an elementary proof that the operator R can entangle quantum states. If 
R is chosen so that ab ^ cd, then the state R(\ip) <8> |^}), with \ip) = |0) + |1), is 
entangled as a quantum state since 

|0) = R(\ij>) ® |^» = a|00) + c|10) + d|01) + 6|11) (6) 

for ab ^ cd (see Fig. 3) is an entangled state. 



2.3 Yang— Baxterization and Hamiltonian 

In this paper, the BGR 6-matrix [""] and the QYBE solution i?-matrix |14[ I15| 
I17| are n 2 x n 2 matrices acting on V V where V is an n-dimensional vector 

space. As b and R act on the tensor product Vi (8> V^+i, we denote them by bi and 

Ri, respectively. 

The BGR 6-matrix has to satisfy the braid relation 

bib i+1 bi = b i+1 bibi +1 , (7) 

while the QYBE has the form 

Ri(x) R i+ i(xy) Ri{y) = R i+ i(y) Ri(xy) R i+ i(x) (8) 

with the asymptotic condition 

R{x = 0) = b, (9) 

and x called the spectral parameter. From these two equations both b and R(x) 
are fixed up to an overall scalar factor. 

In terms of the permutation operator P specified by P(£ <g> rj) = rj <g> £ and the 
i?(x)-matrix, the solution of the algebraic QYBE reading 

Ri2(x)R 13 (xy)R 23 (y) = R23(y)Ri3(xy)R 12 (x) (10) 

takes the form 

R(x) = PR(x) (11) 

where Rij is an operator acting on the tensor product Vi <g> Vj . But we will mainly 
deal with the unitary 4x4 i?(x)-matrix since the unitary permutation operator 
P can be found in order to obtain the unitary i2(x)-matrix. 

As a solution of the QYBE (jHJ), the i?(rc)-matrix usually depends on a deforma- 
tion parameter q and the spectral parameter x. With two such parameters, there 
exist two approaches to solving the QYBE (jHJ). Taking the limit as q — * 1 leads 
to the classical r-matrix satisfying the classical Yang-Baxter equation. Then q- 
deforming r-matrices as solutions of the classical Yang-Baxter equation regain 
the q-dependence for the P-matrices as solutions for the QYBE (jSJ). The QYBE 
solution via this line of thought has been systematized by a general strategy in 
quantum groups (Hopf algebras) approach [T8 l H9l l2?H 121*1 [22j 123] . 

Taking the limit as x — > leads to the braid relation "7J) from the QYBE 
(jHJ) and the BGR 6-matrix from the i?(x)-matrix. Concerning relations between 
the BGR and x-dependent solutions of the QYBE (JHJ), we either reduce a known 
-R(x)-matrix to the BGR 6-matrix, see |25 ( I26| 127] . or construct the R(x)-matrix 
from a given BGR 6-matrix. Such a construction is called Yang-Baxterization. 

In knot theory, these solutions were first studied by Jones [2H an d Turaev 
|27| for the BGR (J7J) satisfying the Hecke algebra relations and for the BGR (|7|) 
satisfying the Birman-Wenzl algebra relations (corresponding to the Kauffman 
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two- variable polynomial, see |2j). Later the more general cases with a BGR (JJJ 
having three or four unequal eigenvalues were considered |281 1291 13(J| , including 
all known trigonometric solutions to the QYBE JHJ). Also Yang-Baxterization of 
non-standard BGR 6-matrix has been discussed [5T1 1521 1551 1531 1551 I5S1 1571 05] . 

In this paper, via Yang-Baxterization the unitary i?(x)-matrices are derived 
and used to construct Hamiltonians determining the time evolution of quantum 
states. In Appendix A, a pragmatic revisit to Yang-Baxterization has been given. 



3 Unitary ^(x)-matrix: the six-vertex model 

In this section, we deal with Yang-Baxterizations of both the standard represen- 
tation and the non-standard representation of the six-vertex model in detail, as 
an example of deriving the unitary ^(x)-matrix with the spectral parameter x 
from the given BGR 6-matrix. 

Consider a non-standard BGR 6-matrix suitable for constructing the Alexan- 
der polynomial [2] 




1 



\ 







,-1 



(12) 



/ q 

1 

\ o o 

where the deformation parameter q has been assumed to be non- vanishing. Let 
b' the transpose and conjugation of b. The unitarity condition b^b = b^b = 1 
leads to q = dbl. 

It has two distinct eigenvalues: q and — q" 1 . By Yang-Baxterization as in 
Appendix A.l, the BGR 6-matrix corresponds to the ^(x)-matrix 



/ q-q' 



R(x) 



V 



: 
(q - q- l )x 
1 — x i 


The case of x 





1 — x 



-i 



\ 



(13) 



qx — q 

1 will not be considered since 



which satisfies the QYBE (| 
R(l) = (q-q- 1 )l. 

Assume that q, x are complex numbers, i / 1 and q ^= 0, with the complex 



conjugations q, x and the norms 
its conjugate matrix by 



|2 ._ 



qq, \\x\ 



f q 
V 



q l x 











(q- q- x )x 


1-x 





1 — X 


q-q- 1 












xx. The i?( re) -matrix has 
\ 









qx 



(14) 



r 1 J 



The unitarity condition shows 

R(x)l&(x) = R\x)R(x) 



pi 



(15) 
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where p is a normalization factor of the #(x)-matrix. It gives us the following 
equations 

' n —1 n2 n —1 Ii2 

\\q — q x\\ = \\qx — q \\ = p 

(1 - x)(q - cT^x + (q - <? _1 )(1 - x) = 

Ik - <? _1 II 2 \\ x \\ 2 + l|l ~ x \\ 2 = P 
\\q — g _1 || 2 + ||1 — x\\ 2 = p 

Comparing the last two equations, we obtain 

1 



1 - q 2 \\ 2 {l - \\x\ 



(16) 



(17) 



which leads to the first case of q = ±1 and the second case of \\x\\ = 1. With 
\\x\\ = 1 and 
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12 ' 



i^ir n<?ir 

going further to rewrite the above second equation as 

1 



(l-x)(l 



r)(9-9)=0 



(18) 



(19) 



which shows that g has to be real. Analyzing the above first equation gives 

7TT2 + <?)(<? - CK^ - = 



(20) 



which is satisfied for real q. To make a safe judgement, comparing the first 
equation and the last two equations, we have 



(q-q)(q(l-x) + q(l-x))=0, 



(21) 



which is also satisfied for real q. Therefore, / 1 and q real, the 

i?(x)-matrix obtained through Yang-Baxterization is unitary. 
In the first case that g = ±1, x ^ 1, we have 



p= ||l - x \\ 2 ■= (l- x )(l-x), 

and then the corresponding ^(x)-matrix is given by 

1 — x 



R(x) 



q=±l 



1 — X 



b\ 



g=±l- 



(22) 



(23) 



In the second case, taking x = e 2ld , 9 7^ and q = e 7 , the _R(x)-matrix takes the 
form 



R(0) = 2 e 1 ' 



/ sinh( 7 - iO) \ 

e l6, simi7 — zsin6> 

— isin<9 e~ l6l sinh7 

\ sinh( 7 + 16) ) 



(24) 
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Modulo the scalar factor 2e td , the unitary i?(#)-matrix has the normalization 
factor p by 

p = sinh 2 7 + sin 2 9. (25) 

In addition, with the new variable 9, such -R(#)-matrix satisfies the following 
QYBE: 



Rni^i + 92) R 



23 i 



(26) 



Rn{0i) -^23(^1 + #2) ^12(^2) — -^23(^1 

Similarly, we can treat the case of q = — e 7 . 

Let consider the BGR 6-matrix which satisfies the Temperley-Lieb algebra 
and plays the role in constructing the Jones polynomial, see |3j, 



/ q 

1 

1 q-q~ l 

\ 






9 ) 



(27) 



It has two eigenvalues: q and — q . With Yang-Baxterization, the corresponding 
.R(x)-matrix has the form 

/ q-q~ l x \ 



R(x) 



\ 



q l x 








(q - q^)x 
1 — x 






1 — x 

i - q 




-1 







q- x x ) 



(28) 



Similar to the preceding case for the Alexander matrix, the unitarity condition 
([15[) informs us that the spectral parameter x lives on the unit circle ||x|| = 1 and 
the deformation parameter q has to be real. 

We apply the ordinary unitarity requirements (|15|) instead of the unitarity 
conditions R(x^ 1 )R(x) oc 1 often mentioned in literature related to the QYBE 
JHJ). With the latter one, we have the following the normalization factor p by 

p = R{x)R{x~ l ) = (q 2 + q~ 2 - 



(29) 



for the non-standard representation and standard representation of the six- vertex 
model. It is compatible with the normalization factor l|25[). 



4 Unitary ^(iz^-matrix: the eight-vertex model 

In this section, we will present the complete solutions of the BGR (J7J) for the 
non-vanishing eight-vertex model and the corresponding unitary i?(x)-matrices 
via Yang-Baxterization. In terms of non- vanishing Boltzman weights w%, u>2, • • •, 
ws, the BGR 6-matrix of the eight- vertex model assumes the form 



( wi 

w 5 

U>4 

\ w 8 





w 3 




w 7 \ 



w 2 ) 



(30) 



9 



Choosing suitable Boltzman weights leads to unitary solutions of the Yang-Baxter 
equation or the braid relation (J7J). In Appendix A. 2, Yang-Baxterization used 
here is sketched in a practical way. 



4.1 Solutions of the eight-vertex model for ^3 = —W4 (I) 



Setting W5 = w\ 



U>2 



wq gives us w{ 



Appendix A. 2. In the case of W3 7^ 104, we have u>3 
BGR ^-matrix has the form 



b± 



( 



\ 












.^1 

U17 








W-j \ 




Wl J 



( 



1 





w\ and w\ + w^ws = 0, see 
and w\ = ±1^3. The 




1 

Tl 






±1 
1 








1 / 



(31) 



It has two eigenvalues Ai = 1 — i and A2 = 1 + %■ The corresponding R{x)- 
matrix via Yang-Baxterization is obtained to be 

-1 



R±{x) 



b± + xAiA 2 b± 

( l + x 



V -T^l- 





l + x 
+(l-x) 






±(1 - x) 
l + x 




\ 



'(1 — X 




l + x / 



(32) 



Assume the spectral parameter x and the deformation parameter q to be 
complex numbers. The unitarity condition (|15|) leads to the following equations 



|1 + x|| 2 + 
11 + x|| 2 + 



II 



x\ 

2 



||1 + x|| 2 + ||1 — x|[ 
(l-x)(l + x) - (l + x)(l-x) 
- x)(l + x) + q (1 + x)(l - x) 

which specify x real and q living at a unit circle. 

Introducing the new variables of angles 9 and (p as follows 

1 x 



P 
P 

P 





(33) 



COS I 



vT+~? 



sm( 



-up 



vT+x 



T 



we represent the i?-|-(x)-matrix in a new form 

R±(0) = cosO b±(<p) + sin0 (6±) _1 (V) 
in which the BGR &±(c/?)-matrices are given by 



1 








1 

+1 







±1 
1 








1 



(34) 



(35) 



(36) 
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4.2 Solutions of the eight-vertex model for W3 = (II) 



Imposing w§ = Wq and W3 = W4 in the given eight-vertex model. In the first case 
of Hi = 2 - t, set z = (t 2 - 2t + 2)5. The BGR 6 ± -matrices take the form 



b± 



( 2 - t q \ 

1 ±z 

±z 1 

\ f 1 t / 



(37) 



It has two distinct eigenvalues 1 ± z and so via Yang-Baxterization the R(x)- 
matrices are obtained by 



R±(x) 



( 2-t(l-x) 



V q'Hl-x) 





1 + x 
±z(l - x) 








\ 



q(l — x t 
±z(l-x) 
1 + x 
2x + t(l-x) J 



(38) 



For simplicity, assume t real and then z positive real. The unitarity condition 
(|T5|) leads to the following equations, 



\\(2 - t) + tx\\ 2 + qq\\l 



x 



q^{{2 -t)+ tx)(l -x) + q(l -x){t + (2- t)x) 
111 + x|| 2 + z 2 ||1 - x\\ 2 



z(l + x)(l - x) + z(l - x)(l + x) 

t)x\\ 2 



rV -1 IU - x\\ 2 + lit + (2 



p 


p 


p 



(39) 



Analyzing the fourth equation leads to ||x|| 2 = 1, while comparing the first equa- 
tion and the third one gives ||g|| 2 = 1. They survive the remaining equations. 
The normalization factor p is given by 



p = 2(1 + z 2 ) + (1 - z 2 ){x + x) = 4 + (t - 1) 2 (2 



x — X 



(40) 



4.3 Solutions of the eight-vertex model for w 3 = w A (III) 

Following the subsection (4.2), consider the second case of W3 = W4, namely take 
±t. The corresponding BGR 6-matrix has the form 



Hi = t so that z 

ui 5 



b± 



( t q \ 

1 ±t 

±t 1 

V q~ l t J 



(41) 



which has three distinct eigenvalues. It needs two types of Yang-Baxterization. 
The first case is essentially the same as the case of two distinct eigenvalues, while 
the second case is treated in detail in Appendix A. 2. 
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With Yang-Baxterization of the first case, we have 

/ t(l-x) q(l + x) \ 

~ . , _ 1 + x ±t(l-x) 

±[X) ~ ±t(l-x) 1 + x 

\ q-^l+x) t(l-x) J 

The unitarity condition (j!5j) shows 

' ||t|| 2 ||l-x|| 2 + gg||l + x|| 2 = p 

\\t\\ 2 \\l -x\\ 2 + ||1 + x|| 2 = p 

||t|| 2 ||l - x|| 2 + ^" 1 ^ 1 ||1 + x|| 2 = p . 

t(l-x)(l+x)+t(l + x)(l-x) = 

q~H(l - x)(l + x) + qi(l + x)(l - x) = 

The first three equations give \\q\\ 2 = 1. Simplifying the fourth equation gives 



(42) 



(43) 



1 — t A 

|2 -i 1 / — \ 

|X|| =1 + - ; o\ X ~ X ) 



where t\ 



Setting t 



|e^ leads to 



x 



l + q 



1 — i tan 4> (x — x) 



(44) 



and setting x = a + ib gives us 



a 2 + 



tan (/)) 2 = sec 2 < 



so that the normalization factor p has the form 

p = 2(1 + ||t|| 2 ) - itan^ (s - s)(l + \\t\\ 2 ) + (1 

Assume t real, namely </ 
given by 



(45) 
(46) 

:\\ 2 )(x + x). (47) 
0, then ||x|| 2 = 1 and the normalization factor p is 

(48) 



p = t 2 (2 -x-x) + 2 + x + x 
which is non-negative since \\x\\ = 1. 



4.4 Solutions of the eight-vertex model for to 3 = w 4 (IV) 

Following the subsection (4.3), the second case for the BGR 6-matrix (|41|) having 
three distinct eigenvalues gives the following i?(x)-matrix: 

/ t(l + x)g X q{l-x) 9l \ 

(l + x)g 2 ±t(l-x)g 2 

±t(l-x)g 2 (l + x)g 2 

\ q- 1 (l-x)g 1 t{l + x) 9l ) 

( t(l + x) q(l~x) \ 

(l + x)g ±t(l-x)g 

±t(l-x)g (l + x)g 

\ q-\l-x) t(l + x) J 



R±(x) 



(49) 
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where the symbols g\, g<i and g are respectively defined by 

gi = l + t + x(l-t), g 2 = l+t-x(l-t), g=—. (50) 

9i 

The unitarity condition Q15jl delivers us the following equations 

f \\t\\ 2 ||1 +x\\ 2 + \\q\\ 2 \\l -x\\ 2 = p 

tq- l {l + x)(l-x)+iq{l-x){l + x) = 

q~ 1 q~ 1 \\l - x\\ 2 + ||t|| 2 ||l +x\\ 2 = p . (51) 

||s|| 2 ||l + x\\ 2 + ||t|| 2 ||5|| 2 ||l - x\\ 2 = p 

i\\g\\ 2 (l + x)(l -x) + t\\g\\ 2 (l + x)(l -x) = 

Identifying the first equation with the third one requires 1 1 1 1 2 = 1. Identifying 
the second equation with the last one shows 

(t - t)(x - x) = 0. (52) 

Focus on the second equation: as t is real, we have ||x|| = 1; as x is real, we have 
t = -i. 

Comparing the first equation with the fourth one leads to the normalization 
factor p. For real t, it takes 

P = \\92\\ 2 = 2(1 + t 2 ) - (1 - t 2 )(x + x); (53) 

for pure imaginary t, it gives 

/3 =|| 52 || 2 = (l-x) 2 + ||t|| 2 (l + x) 2 . (54) 

4.5 Two types of unitarity conditions 

In order to obtain the time-evolution of quantum entangled state determined by 
the unitary ^(x)-matrix, it is better to consider the ordinary unitarity condition 
(fT5|) instead of the unitarity condition R(x)R(x^ 1 ) oc 1. 
For the first type of the eight- vertex model (|32|) , we have 

i?(x)i?(x" 1 ) = 2(x + x- 1 )l, (55) 
which is incompatible with the normalization factor p by 

P = R(x)R\x) = 2(1 + x 2 )1 (56) 

except x = 1. For the second type of the eight- vertex model (|35|) . we have 

R{x)R{x- 1 ) = 2 (1 + z 2 ) + (1 - z 2 )(x + x- 1 ) (57) 

which is compatible with the normalization factor (|40|) for real t. For the third 
type of the eight- vertex model ()42j) . we have 

RtxjRix- 1 ) = 2(1 + t 2 ) + (1 - t 2 )(x + x- 1 ) (58) 

which is obviously an special example of the normalization factor (|47l) . For the 
fourth type of the eight-vertex model (|49j) . we have 

RixjRix- 1 ) = 2(1 + t 2 ) + (t 2 - \){x + x- 1 ) (59) 

which is the same as the normalization factor ()53j) but not (|54(l . 
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5 The i?(x)-matrices as universal quantum gates 



In this section, we will view unitary ^(x)-matrices as universal quantum gates 
with the help of the Brylinski's theorem ^2]. It is a natural generalization of 
the argument regarding a unitary braiding operator as a universal quantum 
gate. Therefore quantum entanglements not only see topological entanglements 
or topological invariants but also know geometric information or geometric in- 
variants hidden in unitary solutions of the QYBE (JSJ) . 
A pure state \ip) denoted by 



1 

i,j=0 



\ l 3) 



® li) 



(60) 



is entangled if aooaii ^ aoiaio- The entanglement of is equivalent to the 
statement that is not the tensor product of two one-qubit states [3D] . 
The .R-matrix having the form 



R 



acts on the tensor product \i) £g> \j) via the formula 



/ R 00 


R° °i 


RT0 




pOl 

"oo 


601 
^01 


601 
K 10 


601 
R n 


510 
-"00 


K>10 

■"oi 


610 
"10 






R 01 


R\o 


m ) 



(61) 



fl|<i> = EE^g|w> 



(62) 



where i,j, k, I take either or 1. The Brylinski's theorem says that it is a uni- 
versal quantum gate when it is a quantum entangling operator which transforms 
the tensor product \tpt P ) into an entangling state R\tptp) given by 

i l i 

R\^ P t) = E E Rihi\ kl ) = E b ki\ kl ) ( 63 ) 

i,j=0k,l=0 



where the coefficients a,j satisfy aooan 
defined by 



k,l=0 

aoiaio and the coefficients b^i are 



>kl 



Efl kl n 

i,j=0 



satisfying b 00 bn b 01 b w . 

Introduce the four dimensional vectors a, b and f. 



ij 



as 



/ aoo \ 
aoi 
aio 
V a n J 



( b 00 \ 
boi 
ho 
V hi ) 



( R 00 \ 

^01 



(64) 



(65) 
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where the upper index T denotes the transpose of the vector so that 



R a = b, 
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(66) 



For convenience of the following discussion, consider the 2x2 matrix A instead 
of the vector a and the 2x2 matrix B instead of the vector b as follows 



A 



aio an 



B 



boo box 
ho b n 



(67) 



The criteria of quantum entanglement for the quantum state determined by 
the matrix A {B) is that the determinant of the matrix A {B) is not zero, namely, 



Det(A) = aooan - a iai ^ 0, Det(B) = b ob n - & i^io / 



(68) 



where the determinant Det(A) (or Det(B)) can be called as the concurrence of 
the corresponding quantum state, see @U1^D- Here in order to judge whether 
the unitary i?(x)-matrix is a universal quantum gate with the Brylinski's theorem 
|13j . we choose Det{A) = (so that the initial state is unentangled) . 



5.1 The case for solutions of the six- vertex models 

For the non-standard representation (|24j) of the six-vertex model, the vector b is 
obtained to be 



( boo \ 
boi 
bio 
\ bn J 



( sinh(7 - i6) a 00 

e ld sinh 7 aoi — i sin 9 aio 
—i sin 6 aoi + e~ te sinh 7 aio 
\ sinh(7 + 16) an 



(69) 



which gives 
boib w 



(sinh 2 7 + sin 2 9) aooan, 



(sinh 2 7 — sin 2 0) aoiaio — isin0sinh7(aoie* 61 + a 2 e 
So the criteria of quantum entanglement has the form 

Det(B) = sin6>[2a oaiisini9 + i(al ie ie + a 2 e _ie ) sinh 7] / 



(70) 



(71) 



Consider the case of the spectral parameter 1 / 1, namely s'm6 7^ 0. The 
choice: aoo = aio = 0, aoi / and 7 / satisfies Det(B) / 0. When 7 = 0, 
namely the deformation parameter q = 1, the criteria of quantum entanglement 
requires aooan 7^ 0. Hence the unitary i?(#)-matrix for the non-standard repre- 
sentation of the six- vertex model is a universal quantum gate except for x = 1. 

For the standard representation (|28j) of the six-vertex model, the determinant 
of the matrix B is a difference between two terms 600^11 an d ^oi&io given by 

^00^11 = sinh(7 — i0) sinh(7 — i6) aoo«ii, 

froi^io = (sinh 2 7 — sin 2 6) aoiaio — isin^sinh7(aoie i6 ' + afge -16 ) (72) 
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The choice: aoo = aio = 0, aoi 7^ and 7 7^ satisfies Det(B) / 0. But the case 
of 7 = leads to Det(B) = 0. Hence the unitary .R(#)-matrix for the standard 
representation of the six-vertex model is a universal quantum gate except for 
x = 1 or q = 1. 

To distinguish the non-standard representation from the standard represen- 
tation in a clear way, the ^(x)-matrices of q = 1 are given respectively by 



Rr 



( 1 \ 

10 
10 

V -1 / 



standard 



( 1 \ 

10 

10 

\ 1 / 



(73) 



in which the non-standard one for constructing the Alexander polynomial is a uni- 
versal quantum gate and the standard one for constructing the Jones polynomial 
is not a universal quantum gate [2]. 



5.2 The case for solutions of the eight-vertex model (I) 

For simplicity of analyzing the criteria of quantum entanglement, we introduce 
the new variable u, v instead of the spectral parameter x, y so that 



1 



1-2/ 



1 — xy u + v 



1 + x ' 1 + y ' 1 + xy 1 + uv 

which suggest the following Yang-Baxter-like equation 



£12 («) £23 ( 



U + V 

1 + uv 



R12M 



1 + uv 



(74) 



(75) 



Specifying a; = e % ' 
forms 



and y = e 1 ® 2 , the parameters in the above equation have the 



u 



-ztan#i, v = —i tan 9 2 , 



u + v 
1 + uv 



-itan(0i + 2 ). 



(76) 



The unitary i?(u)-matrix being a universal quantum gate suggests that the 
unitary ^(x)-matrix is a universal quantum gate since the determinant Det(B) 
are not vanishing in both cases. In terms of the new variable u, the unitary 
R± (it)-matrix for the first type of solution Q32[) of the eight-vertex model has the 
form 

/ 1 
1 





R±(u) 



V 



-q 1 u 









±u 
1 




qu \ 



1 / 



(77) 



The corresponding matrix is given by 



b ± b ± 
°oo "01 



'10 



a o + quan 
=Fu a i + aio 



a i ± u aio 
-q^u a 00 + an 



(78) 
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It has the determinant 

DetiB*) = u(q a 2 u - q~ l a 2 00 ± a 2 01 T a 2 w ) (79) 

which is not zero for the case of it / (i ^ 1), aoo = aoi = and af / aa 2 x . 
Hence the unitary ^-|-(tt)-matrix is a universal quantum gate except u = 0, that 
is to say the unitary i?±(x)-matrix is a universal quantum gate except x = 1. 
In terms of the coefficients a, b, c, d given by 

a oo = ac i °oi = a d, ciio = be, an = bd, (80) 

instead of ciij , the criteria of quantum entanglement leads to 

(d 2 Tq~ 1 c 2 ){qb 2 ±a 2 )^0. (81) 

Consider which type of quantum states violate the criteria of quantum entangle- 
ments. For the fi + -matrix, we have 

d 2 = q- 1 c 2 , or a 2 = -qb 2 , (82) 

while for the R- -matrix, we have 

d 2 = -q~ 1 c 2 , or a 2 =qb 2 . (83) 

Therefore, although a universal quantum gate is identified as an entangling oper- 
ator, it is also able to transform a specified unentangled state (a tensor product 
of quantum states) into another unentangled state. 



5.3 The case for solutions of the eight-vertex model (II) 

In terms of the new variable u, the corresponding R± (u)-matrix has the form 

\ 



R±(u) 



( l + (l-t)u 



V <T X u 




1 







±zu 

1 





qu 



1 + (t - l)u ) 



so that the vector b^ are given by 



V bfi ) 



( (1 + (1 - t)u)a Q0 + qu an \ 
aoi ± zuaiQ 
±zu aoi + oio 
V q~ x u a o + (! + (*- l)u)an j 



which leads to 



z 2 )u 2 )a oan 



(l + (2 

+u(q~ 1 (l + (1 - t)u)a 2 00 + q(l + (t - l)u)a 2 n ), 
(1 + z 2 u 2 ) a iaio ± uz(a^ 1 + a\ 



4) 



(84) 



(85) 



(86) 



To satisfy Det(B) ^ 0, choose u ^ 0, aoo = a oi = a n = but aio / 0. Hence 
the unitary .R-|-(x)-matrix ()38|) is a universal quantum gate for real t and 1. 
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5.4 The case for solutions of the eight-vertex model (III) 

With the new variable u, the R±(u)-matrix is given by 



R±(u) 



( tu q \ 

1 ±tu 

±tu 1 

\ q- 1 tu ) 



(87) 



which says 



(1 + t 2 u 2 ) aooan + ui(ga 2 1 + q 1 d 2 )Q ), 
(1 + t 2 u 2 ) aooan ± ut(oo! + a? ). 



(88) 



So the unitary ^(x)-matrix (fi2*|) is a universal quantum gate except tu = 0. 

The criteria of quantum entanglement in terms of the coefficients (|8fl|) has the 
form 

(^T^q-^qTd^^O. 
For the unitary R± (a;)-matrix, the quantum state specified by 



or 



d 2 



±c 2 q 



(89) 
(90) 



is not entangling even under the action of the unitary -R±(x)-matrix. 

5.5 The case for solutions of the eight-vertex model (IV) 

With the new variable u, the -R±(u)-matrix has the form 



R±(u) 



( t{l + tu) qu{l + tu) \ 

u + t ±tu(u + t) 

±tu(u + t) u + t 

\ q~ l u(l + tu) t(l + tu) J 



(91) 



which gives us 



(t 2 + u 2 )(l + tu) 2 a 00 an + iti(l + tu) 2 (q 1 al + g^n 



(1 + t 2 u 2 ){u + t) z aooan ± ut(u + i)^( a oi + a 



10; 



(92) 



It can be observed that the unitary R± (x)-matrix (f%9"|) is a universal quantum 
gate for ut 7^ 0. 



6 The constructions of the Hamiltonian 

In this section, we present a method of constructing the Hamiltonian from the 
unitary i?(x)-matrix (i?(#)-matrix) for the six-vertex and eight-vertex model. 
The comments on our construction are given in the last subsection. 
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The wave function ip(x) is specified by the unitary i?(x)-matrix, with ip(x) = 
R(x)ip, the pure state ip independent of the time (or the spectral parameter x). 
Hence we obtain the Shrodinger equation corresponding to the time evolution of 
ip(x) controlled by the unitary ^(x)-matrix, 

= H(x)iJ>(x), H(x) = i^^R- 1 (x). (93) 

Here the unitary i?(x)-matrix has to be the form containing the normalization 
factor, namely p^^R{x) satisfying (p^zR)^ 1 = (p~zR)^. The "time-dependent" 
Hamiltonian H{x) is obtained to be 

H(x) = l d -^^{p-lR)-Hx) = i ^p-^t( 4 (94) 
It is expanded into the form 

H(x) = i(d x ]np-*(x) + d x R(x)R' 1 (x)) 
. w ldp(x) dR(x) ~ u 
= lp (*>("2 to + & * t(x)) (95) 

where the formula in terms of w(x) is applied in the following since the calcula- 
tion of R)(x) is easier than that of .R -1 (x). 

In this paper, the spectral parameter x often takes ||x|| = 1, namely x = e td 
so the Hamiltonian H (9) has another form 

m „^ m -im + °m s > m (96) 

which leads to the following Schrodinger equation 

6.1 The case for solutions of the six- vertex models 

For the non-standard representation ()24j) of the six- vertex model, the Hamiltonian 
H(9) is constructed as follows 



H{9) = ^ 



/ coth 7 \ 

— sinh 7 1 

1 sinh 7 

\ -coth 7 / 



(98) 



where the following formulas are used in calculation 



sinh(a + b) = sinh a coth 6 + sinh 6 cotha, sinh ia = i sin a, 
coth(a + b) = coth a coth b + sinh a sinh b, cothia = cosa. (99) 
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Here the Pauli matrices a x , a y and a z are set up as usual 



cr 3 



1 

1 



a 



v 



-i 

1 



a. 



1 
-1 



In addition, the new matrices a± are introduced by a± = \{a x ± icr y ): 

1\ _ / 

ooj' a - ~ [ 1 



OA 



(100) 



(101) 



To represent the Hamiltonian (|98jl in terms of the Pauli matrices, the formulas 
like 



1 



'1 <8 a z + a z 



1) 



( 1 \ 





\ o o o -l / 

are often used. 

The Hamiltonian (|98|) has another form by the Pauli matrices 
sinh 7 



/ \ 

10 



V o o o o J 



(102) 



H{9) 



2p 



■[coth7(l a z + a z ® 1) + sinh 7(1 ® cr a - ^ <g) 1) 



+ (cr + (T_ + (7_ (8) <T+)] 
sinh 7 



2p 



[coth 7(1 (8) (T 2 + <J Z (8) 1) + (cTje (8) CT^ + CJ y (8) CTj, 



+ sinh 7(1 <g> a z - a z % 1)] (103) 

where only the normalization factor p depends on the time variable by 

p = sin 2 6 + sinh 2 7. (104) 

For the standard representation (|28(1 of the six- vertex model, the Hamiltonian 
is obtained to be 



Hie) 



sinh 7 
P 



( coth 7 \ 

— sinh 7 1 

1 sinh 7 

\ coth 7 / 



(105) 



In terms of the Pauli- matrices, it has the form 
sinh 7 



Bid) 



2p 

+ (^+ 
sinh 7 



[coth 7(1 + a z ® a z ) + sinh 7(1 8> <J Z — cr z 
SXT-+CT-® cr+)] 
[coth7l + sinh7(l ® o z — o z ® 1) 



(106) 



+(a x (g) a x + a y (g)(Ty + coth 7 a z (g> a z )] 

which is different from the Hamiltonian (|103j) for the non-standard representation 
of the six-vertex model. 
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6.2 The case for solutions of the eight-vertex model (I) 

With the unitary solution R{x) (|32|) of the QYBE (jHJ), we construct the time- 
independent Hamiltonian H± having the form 



ox 



Interestingly, we have the "time-dependent" Hamiltonian H±(x) by 





( 















i 










Tl 







2 







±1 












V 













/ 



(107) 



H±{x) 



l + x 2 ± 



(108) 



which derives the above Hamiltonian H± at x = 1. When x is real, the Hamilto- 
nian H±(x) is a Hermitian operator. 

For simplicity, we study the unitary i?(#)-matrix (|35j) to decide the unitary 
evolution of quantum states. After some algebra, the Hamiltonian H± is obtained 
to be 



ldp(e) , 3fl ± (S) B 



ff ±w = i„-i( e0 (-- ir + 

which is independent of the time variable 9. 

In terms of the Pauli matrices a x , a y and a z and a±, the Hamiltonian (|109|) 
has the form 



± 2 



cr + + e' v cr_ (gi a_ =F <x+ ® a- ± a. 



(110) 



Introducing the two-dimensional vector a and two unit directional vector Hi 
and n,2 in icy-plane: 



, , ^ , ir + ip . vr + y? _ </9 . p 

cr = {a x ,a y ); n x = (cos — - — ,sm — - — J, n 2 = (cos — , sin —J 



the projections of the vector a into n\ and n2 are given by 
a 



T II2 



a • ni = a+e'i^ + a.e^ +7T \ 



(t„„ — a ■ n2 = a + e ^ + a -e^ . 
The Hamiltonian (jllOj) can be recast to 

if. 



1 



-a ni <g) cj n2 , 



1 



(111) 



(112) 



(113) 



Consider the time-evolution operator U±{9) determined by the Hamiltonian 
H±, for example, U+{9) given by 



U+(9) 



e 2 



(o- ni (g>o-, l2 )0 



COS ■ 



i sin - a ni 



(7 



n 2 - 



(114) 



21 



6.3 The case for solutions of the eight-vertex model (II) 

The unitary R±( x ) -matrix (|38[) requires x = e td , ||g|| 2 



1 and real t. The 



normalization factor p and its derivative have the forms 



p = 4 + 4(t - l) 2 sin 2 



(9/9 



2(t- l) 2 sin<9. 



2' <9# 

After some calculation, the Hamiltonian H±{6) is obtained to be 



(115) 



H±(9) 



ldp{9) dR±{9) 



Rim 



2 89 89 

( 1 - 1 q \ 
±z 







dbz 







V 



(116) 



t-1 J 



In terms of the Pauli matrices, it is shown up as follows 

H±(0) = -\ + {l-t)p- l (\®a z + a z ®t) 

2 i 
H — \qa + <g) a + + q cr_ <g> <r_ ± z(cr + cr_ + <r_ er + ) . (117) 
P 

Consider two special cases of making the normalization factor p independent 
of the time variable 9. Take 9 = 0, namely, choose the Hamiltonian H±(9) defined 
as 

H±(9) = i-^(p- 1 iR ± (9))\ e=0 (118) 
which leads to the time-independent Hamiltonian 



H±\e= 



1 



[qa+ ® a + + q 1 cr_ <g> cr_ ± z(a+ cr_ + <r_ (g) cr + )] 



-1 + ® + ff* ® 1). 



(119) 



Take the parameter t = 1, the time- independent Hamiltonian //-t is given by 



i/± = -[— 1 + <7<r + <8> + + g x O- <r_ ± (cr + <g) <r_ + <r_ (8) 0+)]. 
With the two unit directional vector ni and n2 in xy-plane: 



ni = (cos-, sin-), n 2 = (cos — — ,sin— — ), 



the Hamiltonians ()12U|) have the following forms 

H + = ~(-l + a ni ®a ni ), H- = --(! + an* 
Consider the unitary time-evolution operator U±(9), for example, U+{9) given by 



(120) 
(121) 
(122) 



U+(9) 



-iH 4 



e 2 (cos - - i sin - a ni a ni ] 



(123) 
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6.4 The case for solutions of the eight-vertex model (III) 

Choose real t, x = e lS and \\q\\ = 1 for the unitary /2-|-(x)-matrix (|42[) . The 
normalization factor p and its derivative are given by 



p = 4 + 4(t 2 -l)sin 2 -, 



dp 



2' 86 

The corresponding Hamiltonian H±{6) is found to be 



2(i 2 - l)sin0. 



(124) 



1. 



-1 + 





( 








q \ 


2t 








±1 





P 





±1 










V q- 1 








o / 



(125) 



H±(9) 

which leads to the Hamiltonian represented by the Pauli matrices 

1 2t 

H±(6) = --1 + —[qa + (g> a + + q~ X (J- ®(J-± (a + ® cr_ + <r_ ® a+)]. (126) 

2 p 

In the case of = 0, the time-independent Hamiltonian is given by 
H±\ s=0 = --1 + ^[gcr+ (8) cr + + g _1 cr_ (8> cr_ ± (cr + (8 cr_ + <t_ (8 cr + )]. (127) 

In the other case of t = 1, the other time- independent Hamiltonian H± has the 
form the same as (jl2U|) . 

6.5 The case for solutions of the eight-vertex model (IV) 

Consider x = e lS , \\q\\ = 1 and real t for the unitary A^(a;)-matrix (|49|). The 
factor ||<?i|| 2 , 1152 1| 2 and the normalization factor p take the forms 



1 91 1 



2(1 + t 2 ) + 2(1 - t 2 ) cos (9, ||g 2 II 2 = 2(1 + t 2 ) - 2(1 - t 2 ) cos i 

■ 2 ■ I „ || 2 An , *2\2 .i/i .2\2 „ 2 , 



|2 

P = = 4(1 + 4(1 -t*) 2 cos* 0. 

Through some calculation, the Hamiltonian H±{6) is obtained to be 



(128) 



2t 
P 



/ IM| 2 










llfflll 

M\92\ 



H ± (0) 

P \ q- l \\9if 

which gives the following Hamiltonian as 





±11521 

\\gif 




Of 1 1 

H± (6) = -1 + -[-(||5i|| 2 + IM| 2 )1 ® 1 + -(||ff 2 || 2 
p 2 2 



Ibill 2 \ 



Ml 2 / 



bill 2 )^ ® ^ 



(129) 



+ ||gi|| 2 (gcr + ®cr + +g V_ (8) cr_) ± |M| 2 (o-+ ® o - - + (8) a+)]. (130) 
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The time-independent Hamiltonian obtained by taking 8 = is given by 

H±\ e=0 = -1 + ^-[(t 2 + 1)1 (8) 1 + (t 2 - l)a z <g> a z 
+2(qa + <g> a + + q~ X G- <g) <?_) ± 2t 2 (a + <g> cr_ + cr_ cr+)], (131) 

while the other time-independent Hamiltonian H± obtained by taking t = 1 is 
the same as (jl2()|) . 

6.6 Comments on our constructions of Hamiltonian 

The Schrodinger equation is a differential equation. It's solution represents the 
evolution of the initial state (input). Here before treating the Schrodinger equa- 
tion, the time evolution of the state (the unitary braiding operator) is known as a 
discrete evolution and is determined by the unitary i?(x)-matrix. Our problem is 
to find out which type of the Schrodinger equation has evolutionary solutions the 
same as the evolutions given by the unitary i?(x)-matrices. In our case, therefore, 
the Schrodinger evolution is recognized as the unitary ^(x)-matrix. 

However, the Schrodinger equation says more than just the unitary R{x)- 
matrix. It leads us to study the physics behind it, and gives us new unitary 
solutions which are not necessarily solutions to the QYBE. The construction of 
a Hamiltonian from the unitary ^(x)-matrix shows that the QYBE is a physical 
subject. In history, the Schrodinger equation is the starting point and how to 
solve the Schrodinger equation is a central topic. Now we try to recover (relate) 
physics from (to) the unitary ^(x)-matrix. 

In our case, the construction of the Hamiltonian in terms of the unitary R(x)- 
matrix follows a traditional approach in literature on the QYBE. There are two 
ways of choosing the time- variable. Each choice is determined by a corresponding 
purpose. We explain the spectral parameter x (6) as the time variable. That is 
to say that we choose the time-evolution of quantum state (the unitary braiding 
operator) as the unitary ^(x)-matrix. This is a natural and necessary choice in 
the sense of regarding unitary i?-matrices as universal quantum gates. We want 
to see the evolution of unitary braiding operator or quantum state. 

Indeed in the literature, the ordinary time in space-time as the time variable, 
the spectral parameter can be explained as the momentum, but the time-evolution 
of quantum state determined by the Hamiltonian can not be identified with the 
unitary i?(x)-matrix. The physics in the case is well-known such as XXX model 
(and so on). 

Finally it seems that the Schrodinger equation in our construction does not 
have space variables like the ordinary Schrodinger equation has. We could explain 
that we have physics on lattices of space: discrete physics. It seems that physics 
in our case is close to that on spin chains like XXX model (and so on). 
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7 The CNOT gates via the i?-matrix 



The gate G is universal for quantum computation (or just universal) if G together 
with local unitary transformations (unitary transformations from V to V) gener- 
ates all unitary transformations of the complex vector space of dimension 2 n to 
itself. It is well-known [l] that the CNOT gate is a universal gate. 

In |10| . Kauffman and Lomonaco prove the following result. 

Theorem 1. Let 



R 



( l/y/2 l/y/2 \ 

l/y/2 -1/V2 

1/V2 1/V2 

V -l/y/2 l/y/2 J 



(132) 



be the above unitary solution to the braid relation (J7J). Then R is a universal 
gate. The proof below (repeated from ^Hj) gives a specific expression for the 
CNOT gate in terms of R. 

Proof. This result follows at once from the Brylinksis' theorem ^3], since R is 
highly entangling. For a direct computational proof, it suffices to show that the 
CNOT gate can be generated from R and local unitary transformations. Let 

l/y/2 1/V2 \ (-l/y/2 l/y/2 

l/y/2 -l/y/2 J' 1 \ i/y/2 ijy/2 

l/y/2 i/y/2 \ ( 1 

l/y/2 -i/y/2 )' \ i 



o 



(133) 



Let M = a <g> /3 and N 



This completes the proof. 



-7 (8> 5. Then it is straightforward to verify that 
CNOT = M-R- N. 



□ 



We now show how Yang-Baxterization illuminates the structure of this Bell 
basis transformation. We discuss physics related to the time-evolution of the 
universal quantum gate determined by the unitary .R(#)-matrix ()35|) . The braid 
group representation 6±((/?)-matrix (|36|) yields the Bell states with the phase 
factor e lip , 



b±(<p) 



( |oo) \ 
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|00) -e^|ll) 
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|01) 


1 




|01> =F |10> 




|10) 






±|01) + |10) 




V in) ) 




\ 


e-^|00) + |11) 


/ 



(134) 



which shows that ip = leads to the Bell states, the maximum of entangled states, 
-±=(|00)±|11», -^=(|10)±|01)). (135) 
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In terms of the Hamiltonian H± (|109[) . the P±(#)-matrix (j35[) has the form 



R±{0) = cos 



7T 



7T 



+ 2i sin(- -6)H± 



(136) 



which can be also used to construct the CNOT gate with additional single qubit 
transformations, examples see ^0]- The unitary P-matrix (|132j) is realized by 



R — R-{9)\e =ip =Q 
In addition, with the unitary P(#)-matrix (j35j) . we have 



(137) 



R±{0) 



( |00) \ / cos(f -0)|OO) -e^sin(f - 9)\ll) \ 
|01) = cos(| -0)|Ol)Tsin(f -0)|1O) 
|10) cos(i-0)|lO)±sin(f-0)|Ol» 

\ |H) / \ cos(f-0)|ll) + e -^ S in(! -0)|OO> / 



(138) 



Hence with the concept of the Bloch vectors on the Bloch sphere £Q , the variables 
9 and ip realize their geometric meanings and the construction of the CNOT gate 
becomes clear. 

To obtain other two-qubit quantum gates, for instance the CNOT gate, we 
have to apply single qubit unitary transformations A,B,C,D which can be pos- 
sibly found in the Bloch sphere [Q by SO(3) rotations, namely, 



(A ® B)U±(0)(C ®D) = P 1 ®1 + P l ®a x = CNOT 



(139) 



which yields the CNOT gate and where the states | f ) and | j) are the eigenvectors 
of a z , a z \ j) = | T),0z| I) = — I I) an d the projection operators Pf and Pj_ have 
the forms 



1 



TXT I, 



Pi 



I 



Define the SO(3) rotation around the n-axis by 



D n (0) 

where a = (a x , a y ,a z ). For examples: 



e 2 



{3-n) 



D z (-^)=e^% D x 



■7T> 

T 



satisfy 



C^)P',(-|K 1 P» z (|)P» s (--)=a a 



^(-|)o-n 2 ^(|) = ^ 



(140) 

(141) 
(142) 
(143) 



Consider the time-evolution operator U+(0) (|114[) . Choosing suitable single 
qubit transformations, we obtain 



. 7T 



(D X (-)D Z (-^-) ® d s (-%))U+(P)(D z {%)D x (-~) ® £>*(£)) 



e 2 



(144) 
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which has another form 

e -4(<T,ia<T«)0 _ p T e -5 CT - e + p i o e 3»-9. (145) 

Set 9 = | . To construct the CNOT gate, we need additional single qubit trans- 
formations 

(S <g> e i i ff ")e- < T( <T *® '") = CNOT (146) 

in which the phase gate 5 has the form 5 = Pf — i Pj_ , see (|133[1 . 

Consider the time-evolution operator £/_(# = — 5, <£> = 0), namely the unitary 
P-matrix fT32}) given by e 1 ^*®^ which is transformed into e *f ( (T *® <r *) by 

(A,(-|) ® jDz (-|)) e ^K^)( J D y (|) g, = e ^(-^). (147) 

So we obtain another proof for Theorem 1 jl()j . 

8 Concluding remarks 

Motivated by the observation that there are certain natural similarities between 
quantum entanglements and topological entanglements, we derive the unitary so- 
lutions of the QYBE or the unitary R{ x ) -matrices via Yang-Baxterization and 
construct the related Hamiltonians for the standard and non-standard represen- 
tations of the six-vertex model and the complete solutions of the non- vanishing 
eight-vertex model. With the Brylinksis' theorem .13., the unitary P(j;)-matrix 
is also a universal quantum gate except very special cases. 

The remark has to be made on the classification of the unitary P(x)-matrices 
we obtained. We classify them according to Yang-Baxterization. We don't try to 
classify them with other approaches since we can't obtain the complete solutions 
of the QYBE (jHJ) only via Yang-Baxterization. But the unitary solutions of the 
Yang-Baxter equation (the braid relation) have been classified F° r example, 
the first type of solution for the eight-vertex model in our case belongs the fourth 
family specified by But for the eight- vertex model, we have three types 

solutions of the BGR Q but have four types of the unitary P(x)-matrices. 

It is worthwhile arguing again that quantum entanglements for quantum infor- 
mation processing are related to not only topological entanglements but also geo- 
metric invariants. The previous one leads us to viewing the unitary braiding op- 
erator as a quantum entanglement operator and a universal quantum gate, while 
the latter one suggests us to regard unitary P(x)-matrices as universal quantum 
gates. In addition, it is important to mention another view of topological issues 
for quantum computing in terms of anyonic models, see |42 [ I43| IH | I45 | HH j 147] . 
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A A practical revisit to Yang— Baxterization 

Yang-Baxterization [21] is a prescription deriving solutions to the QYBE 
from the BGR 0. We consider a BGR fe-matrix with two distinct non- vanishing 
eigenvalues Ai and A2 taking the form 

b = Ai Pi + A 2 P 2 (148) 

where P\ and P 2 are the projection matrices satisfying 

Pi + P 2 = 1, Pi 2 = Pi, P 2 2 = P 2 , Pi P 2 = 0. (149) 

With the help of Yang-Baxterization, the corresponding P(x)-matrix modulo an 
overall scalar factor |281 129j has the form 

R(x) = (Ai + A 2 x) Pi + (A 2 + Ai x) P 2 

= 6 + A1A2X6" 1 (150) 

where the inverse matrix b^ 1 is given by 

& -1 = T P i + T p 2- (151) 
Ai A2 

In the case that the BGR 6-matrix has three distinct non- vanishing eigenvalues 
Ai, A 2 and A3, Yang-Baxterization |28[ I29j leads to the following formula 

R(x) = AiA 3 a;(x - 1)6 _1 + (Ai + A 2 + A 3 + A1A3A2 l )xl - (x - 1)6. (152) 

Changing the ordering of three different eigenvalues in principle gives us different 
solutions of the QYBE (jHJ). However, this formula is symmetric with respect to 
interchanging Ai and A3, so usually we obtain three types of the P(x)-matrices. 
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With two formulas (j!50[) and (|152[) . the i?(x)-matrix is proportional to the 
unit matrix 1 at x = 1, namely, 



R(x 



1) oc 1 

-l 



(153) 



except the case of Ai + A2 + A3 + AiAaA^ = 0. Hence for simplicity, in the 
following we will not pay special attention to the case of x = 1. 

One remark has to be made. The formula (j!5(Jj) has been proved to satisfy the 
QYBE (JHJ) , however the formula (|152j) does not have the general proof verifying 
that it is the solution of the QYBE (JHJ). Once the ^(x)-matrix is obtained via 
Yang-Baxterization, it would be safest to check whether it truly satisfies the 
QYBE (JHJ). 

A.l Yang— Baxterization of the six- vertex model 

Consider a non-standard BGR 6-matrix suitable for constructing the Alexander 
polynomial [2] 



\ 



,-1 



(154) 



/ q 

1 

1 q - q' 1 

\ -q 

where the deformation parameter q has been assumed to be non- vanishing. It 
has two distinct eigenvalues: q and — q~ l . In terms of the projection matrices 
Pi(g) andP 2 (g): 
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1 ) 



(155) 



the BGR 6-matrix is also given by 

b = qP 1 (q)-q- 1 P 2 (q). 



(156) 



With the help of Yang-Baxterization, the BGR 6-matrix corresponds to the fol- 
lowing i?( re) -matrix satisfying the QYBE (JHJ), 

(q - q- 1 x)P 1 (q) + (-q' 1 + qx)P 2 (q) 
b-xb- 1 . (157) 



R(x) 



A. 2 Yang— Baxterization of the eight-vertex model 

The eight-vertex model assumes a general form 
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w 7 \ 
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w 2 J 



(158) 
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in terms of non- vanishing Boltzman weights uij, i = 1, • • • 8. We will present the 
complete solutions of the BGR ((7J) for the non-vanishing eight-vertex model and 
construct the corresponding fi(x)-matrices via Yang-Baxterization. 
Introduce by 



hup 



Xlu 
HP 



b pv b Xp b UK ) 



(159) 



where indices take spin up or spin down, namely ± sign or spins ±|. With 
the above ansatz, we have b^i = for a + (3 ^ k + to mod 2. In the case of 
a + /3 + 7^ A + u> + k mod 2, the braid relation ® is satisfied automatically, 
= ^- Hence we have to treat thirty-two equations of the Boltzman weights. 
We solve all possible equations obtained by substituting the above ansatz 
of the eight-vertex model into the braid relation (|7j). We observe the equation 
(u>5 — wq)wjW8 = showing 105 = wq. With the equations (103 — Wi)(wi — w^wg = 
and (1V3 — W4)(w2 — w^)wj = 0, we have to choose either W3 = W4, w§ = w§ or 
w 3 / W4, w 5 = wi = w 2 = Re- 



setting 1V5 = w\ = W2 = wq gives us w\ 



the case of w% ^ W4, we have 103 
have the forms 



-W4 and w\ 



= w\ and w\ + w-jwg = 0. In 
±W3. The BGR 6±-matrices 



6± 
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W\ u>7 \ 
ioi ±t«i 
w\ 



=fri 



U17 



\ 



1 q \ 

1 ±1 

=Fl 1 

1/ 



(160) 



It has two eigenvalues Ai = 1 — i, A2 = 1 + i. 
via Yang-Baxterization are obtained to be 



The corresponding ^(x)-matrices 



R±(x) = b± + 2xb± l . 



(161) 



Imposing w§ = wq and W3 = in the given eight-vertex model, the solutions 
of the QYBE © have to satisfy the following three independent equations 29 , 



105 - W 7 W 8 
w 2 ~ w 3 + ^1^5 — 



™1 



0, 

0. 

0. 



(162) 



Since the i?(x)-matrix is modulo an overall scalar factor, we introduce 



q = — , q = — , t = — , 
and rewrite the above equations in an explicit way as 



W3 
w 5 



(163) 
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Therefore we obtain two types of the BGR (|7)l. The first one is given by taking 
= 2 — t and the second one is specified by ^ = t. 

1 

In the first case of — 
take the form 



2 - i, set z = (t 2 - 2t + 2)2. The BGR 6±-matrices 



b± 



( 2-t 

1 

±2 

V q' 1 



g \ 
±z 

1 
t ) 



(165) 



It has two distinct eigenvalues 1 ± z and so via Yang-Baxterization the R{x 
matrices are obtained by 



i?±(:c) = 6± + x(l 

where the matrix entries satisfy W1W2 + w^w^ = w^wq + wtw%. 

In the second case of ^3 = W4, namely take ^ = t so that 
corresponding BGR 6-t-matrix has the form 






Vt 1 
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±t 





g \ 
±t 

1 
t ) 



(166) 



±t. The 



(167) 



It has three eigenvalues 1 + 1, 1 — t and — 1 + 1. To apply the formula (|152j) for 
the case of three distinct eigenvalues, we arrange three eigenvalues Ai, A2 and A3 
in three different orderings as follows 



Ai A2 A3 

the first ordering: 1 + 1 1 — t t — 1 

the second ordering: 1 +t t — 1 1 — t 

the third ordering: 1 — t 1 + t t — 1 



(168) 



For the first ordering and the second one, we have 

A1 + A2 + A3 + A1A3A2 1 =0 (169) 

so that the situation goes back to the case of two distinct eigenvalues, 

R ± ( x ) = -(x-l)(6±± x(l-t 2 )b^) 

oc b± ± x (1 - t 2 ) 6 ± x (170) 

where the plus is for the first ordering and the minus is for the second ordering. 
We only need discuss the case of the first ordering. 

For the BGR 6-matrix (|167j) having three distinct eigenvalues, with the first 
ordering, we have 

(171) 
(172) 



t 1 



S 1 ; 



R±(x) = b± + x(l 
with the third ordering (|168|) . we have 

R±{x) = -{x- l)X 2 (b± - A1A3X bg 1 ) + {\ 2 2 + AiA 3 )xll. 
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